Widefield in vivo imaging system with two fluorescence and two reflectance channels, a single sCMOS detector, and shielded illumination

Abstract. Significance Widefield microscopy of the entire dorsal part of mouse cerebral cortex enables large-scale (“mesoscopic”) imaging of different aspects of neuronal activity with spectrally compatible fluorescent indicators as well as hemodynamics via oxy- and deoxyhemoglobin absorption. Versatile and cost-effective imaging systems are needed for large-scale, color-multiplexed imaging of multiple fluorescent and intrinsic contrasts. Aim We aim to develop a system for mesoscopic imaging of two fluorescent and two reflectance channels. Approach Excitation of red and green fluorescence is achieved through epi-illumination. Hemoglobin absorption imaging is achieved using 525- and 625-nm light-emitting diodes positioned around the objective lens. An aluminum hemisphere placed between objective and cranial window provides diffuse illumination of the brain. Signals are recorded sequentially by a single sCMOS detector. Results We demonstrate the performance of our imaging system by recording large-scale spontaneous and stimulus-evoked neuronal, cholinergic, and hemodynamic activity in awake, head-fixed mice with a curved “crystal skull” window expressing the red calcium indicator jRGECO1a and the green acetylcholine sensor GRABACh3.0. Shielding of illumination light through the aluminum hemisphere enables concurrent recording of pupil diameter changes. Conclusions Our widefield microscope design with a single camera can be used to acquire multiple aspects of brain physiology and is compatible with behavioral readouts of pupil diameter.


Introduction
3][4][5][6][7] Probes available in distinct spectral variants 8 can be spectrally multiplexed for simultaneous imaging of different aspects of neuronal activity. 9Using multiplexed probes for in vivo widefield imaging of head-fixed mice with "crystal skull" cranial windows 10 or transcranial imaging, optionally combined with thinned-skull preparations, 11 provides a view of the entire dorsal part of the cerebral cortex including the primary sensory, motor, and higher association areas. 12,13n the healthy brain, increases in neuronal activity are often accompanied by increases in blood flow, volume, and oxygenation in the active region.The increase in blood flow and volume occurs, in large part, due to dilation of cerebral arterioles responding to vasoactive messengers released by active neurons. 14,15This hemodynamic response increases oxygen supply beyond the actual demand leading to a net increase in blood oxygenation.Both oxy-and deoxyhemoglobin (HbO and HbR, respectively) strongly absorb light in the blue-green spectrum.Therefore, the hemodynamic response results in time-variant modulation of fluorescence signals within that spectral range.The methods developed to correct fluorescent signals for this hemodynamic artifact fall into two categories.For the first correction method, enhanced green fluorescent protein (EGFP)-based indicators, including GCaMP, iGluSnFR, and GRAB indicators, are excited off-peak at their isosbestic point in the 390-to 420-nm range.When excited at these wavelengths, fluorescence emission is independent from the analyte concentration (i.e., ½Ca 2þ , [glutamate], etc.) but affected by hemodynamic changes. 16,17This signal can then be regressed out of the analyte-dependent signal to remove hemodynamic artifacts.Correction methods measuring changes in hemoglobin absorption use either a single reflectance measurement at 520 to 530 nm or estimate changes in HbO and HbR concentrations (Δ½HbO and Δ½HbR) from simultaneous acquisition of light absorption from at least two wavelengths. 18Furthermore, quantitative estimates of Δ½HbO and Δ½HbR are important measurables on their own; they aid in understanding cerebrovascular physiology and interpretation of noninvasive imaging signals generated through, for example, functional magnetic resonance imaging. 19urrent stationary widefield imaging setups typically use light-emitting diodes (LEDs) with sequential illumination, one or two low-noise, high-sensitivity scientific CMOS cameras, and low (smaller than 2×) magnification optics with a large field of view. 20Experimental demands led to various adaptations of this basic blueprint.In systems dedicated for widefield calcium imaging with EGFP-based indicators, blue LEDs (450 to 480 nm) are used for fluorescence excitation, with light delivered through the objective lens or through side-illumination in some cases.0][31] Several groups incorporated a second fluorescence channel with excitation at 560 to 590 nm to enable two-color fluorescence imaging of red fluorescent calcium indicators like jRCaMP1b or jRGECO1a in combination with either EGFP-derived indicators, such as GCaMP6, 32 GRAB ACh , 33 or flavoprotein autofluorescence. 34,35Widefield imaging systems dedicated to the investigation of neurovascular interactions typically combine imaging of green 4,7,30,[36][37][38] or red fluorescent indicators, 34,39,40 or green and red fluorescent indicators 35 with reflectance imaging for Δ½HbO and Δ½HbR estimation at two or three wavelengths, typically in the range of 525 to 565 nm and 590 to 625 nm.Light for reflectance imaging is delivered directly from one or more multicolor light engines, 7,30,40 through a single light guide, 39 through an illumination ring, 22,41,42 or through several optical fibers connected to a single light source 38 to illuminate the sample surface at an angle.Some researchers report the use of linear polarizers in front of the light sources to minimize specular reflection from the glass window. 7,30,42For experiments involving behavioral tasks or visual stimulation, different solutions to block illumination light from reaching the animal's eyes have been reported, for example, through implanted light shielding, 38,43 custom light-shielding between objective and headplate, 3,44 or imaging stage design. 45ere, we describe a versatile widefield imaging system that enables two-color fluorescence imaging in addition to concurrent acquisition of two reflectance channels for quantitative estimation of Δ½HbO and Δ½HbR with a single sCMOS detector.Specular reflections from curved glass cranial windows that interfere with reflectance imaging are minimized using diffuse illumination with 525-and 625-nm light for reflectance imaging by combining a custom-built ring illuminator with an aluminum hemisphere between the objective and the cranial window.In addition, the aluminum hemisphere shields illumination light to prevent it from reaching the animal's eyes and thereby avoids unwanted visual stimulation of the animal.We provide optical design blueprints, a wiring diagram, a timing chart, data acquisition software, and an inventory of all parts.We demonstrate the performance of our new widefield microscope in awake mice expressing the mApple-based calcium indicator jRGECO1a 46 and the EGFP-based acetylcholine (ACh) probe GRAB ACh3.0 . 17

Results
We designed and manufactured a widefield imaging system with a field of view of 10 × 10 mm 2 covering the dorsal cortex of a mouse.We use a single sCMOS detector for sequential recording of two fluorescent indicators (here EGFP-derived GRAB ACh3.0 and mApple-derived jRGECO1a) as well as hemoglobin absorption at 525 and 625 nm at an effective frame rate of 10 Hz or better without changes in the configuration of the emission path.Figure 1 and Fig. S1 in the Supplementary Material provide an overview of the system layout and illustrate the spectral properties of the light sources and filter elements together with HbO and HbR absorption spectra as well as jRGECO1a and GRAB ACh3.0 excitation and emission spectra.For fluorescence excitation, we use epi-illumination through the objective lens.A 470-nm LED with a maximal power of 3.7 W provides light for GRAB ACh3.0 excitation, which is filtered through a 466/40 nm excitation bandpass filter, resulting in a maximal power of 403 mW at the sample surface.For jRGECO1a excitation, we use light from a 565-nm LED with a maximal power of 8.6 W, which is filtered through a 560/14 nm excitation bandpass filter reaching a maximal power of 75 mW at the sample surface.A secondary peak in the 565-nm LED spectrum at ∼425 nm is removed by an additional 450-nm longpass filter in the excitation path [Fig.1(a)].The collimated and filtered light is coupled into the infinity space using a multiband dichroic mirror (488/561 nm, "main dichroic").This, in combination with a multiband "detection filter" (523/610 nm), allows GRAB ACh3.0 and jRGECO1a fluorescence at 500 to 540 nm and 580 to 640 nm, respectively, to reach the sCMOS detector while blocking excitation light at 440 to 480 nm and 550 to 570 nm from reaching the camera [Figs.1(d), 1(e), and 1(g)].Typical powers at the sample surface during image acquisition were 90 to 177 and 68 to 75 mW for 470 and 565 nm, respectively, at illumination times of 3 to 6 ms.We observed no or minimal bleaching over continuous acquisition periods of up to 30 min.The wavelengths for hemoglobin absorption imaging are within the respective emission bands of the green (525 nm) and red (625 nm) fluorophore [Fig.1(f)].The LEDs used for hemoglobin absorption imaging are placed in a custom-designed reflectance illumination ring (RIR) located between the objective lens and the sample surface [Fig.1(b) and Fig. S1(D) in the Supplementary Material].The RIR contains six 525-nm LEDs and six 625-nm LEDs positioned at equal distances with maximal powers of 170 mW and 371 mW per LED, respectively.To enable a simple procedure for estimation of Δ½HbR and Δ½HbO, we placed 10-nm bandpass filters in front of each LED.This leads to a maximal power of 39.5 mW (525 nm) and 51 mW (625 nm) at the sample surface.During image acquisition, typical powers at the sample surface were 10.5 and 4.4 mW for the 525-and 625-nm LEDs, respectively, at 3-ms illumination times.We placed a hollow aluminum hemisphere between the RIR and the sample surface internally reflecting the light of the 525-and 625-nm LEDs.For imaging, the opening of the hemisphere is brought close to the sample surface achieving diffuse illumination and efficiently eliminating most reflections from the curved glass window covering the brain surface.
With the chosen sCMOS camera working in rolling shutter mode, we ensure that the LEDs emit light only when all rows of the sCMOS chip are simultaneously exposed.Therefore, we control the LED drivers with two digital trigger signals combined through a logical AND gate integrated circuit (IC) [Fig.1(c)].The first input to the AND gate IC is generated by the camera and is in "high" state when all rows are being exposed; the second input is generated by the microscope control computer and selects the respective LED wavelength to be switched on [Fig.1(c)].Only when both signals are in high state, the LED power supply receives a trigger The effective frame rate is determined by the length of the four-frame cycle (typically 100 ms).Grey lines next to the camera Fire ALL trigger represent exposure of pixel rows in rolling shutter mode.The Fire ALL trigger is sent when exposure overlaps for all rows.Each LED driver (bottom left) is controlled by the output of a separate AND gate.The Fire ALL trigger is the common input to every AND gate, whereas the individual LED selection triggers are the second input for each AND gate.(d).Top: light-blue and dark-blue shaded areas represent the 470-nm LED spectrum before and after passing the excitation filter (BP466/40, blue dotted line).Bottom: excitation (solid blue line) and emission (solid green line) spectrum of EGFP-derived GRAB ACh3.0 . 17(e) Top: light-orange and dark-orange areas represent the 565-nm LED spectrum before and after passing the excitation filter (BP560/14, orange dotted line).Bottom: excitation (solid yellow line) and emission (solid red line) spectrum of mApple-derived jRGECO1a. 16f) Top: LED spectra are confined by excitation filters (525 nm in light and dark green; 625 nm in light and dark red before and after passing through the respective excitation filters).Bottom: absorption spectra of oxyhemoglobin (HbO, red line) and deoxyhemoglobin (HbR, blue line) from Prahl.47 (g) Spectra of the main dichroic (brown) and detection filter (pink).Light must pass through both filters to reach the camera.
signal to switch the respective LED on.Timing delays of the AND gate IC and the LED driver are <9 ns and <100 μs, respectively, enabling virtually delay-free LED power switches for typical exposure times of 2 to 8 ms per individual LED.
The imaging setup is controlled through a DAQ interface programmed via MATLAB.The main acquisition script generates digital trigger signals controlling the imaging setup based on exposure time, which can be different for every wavelength, acquisition rate, and total acquisition time.Additional triggers are generated to synchronize or control external devices, such as stimulus apparatuses, behavioral cameras, electrophysiology recording equipment, and behavioral control devices.The currents of the individual LEDs are set in MATLAB and sent via a USB-based remote-control interface to the LED drivers.The sCMOS camera is controlled through the manufacturer's software, where parameters, such as field of view and on-camera pixel binning are set.This software also drives image acquisition; however, the camera is set to respond to exposure triggers generated by the DAQ interface, with the (variable) trigger length determining the exposure time.During acquisition, the DAQ interface records feedback signals from camera, microscope control, and external devices allowing synchronization of data streams during off-line postprocessing and analysis.
Although our single-camera approach reduces complexity and cost, we use multiband filters to image two fluorophores present in the sample at the same time.The lack of spectral separation on the emission side can potentially create crosstalk when fluorophores are excited offresonance.Note that both fluorophores undergo dynamic changes in fluorescence intensity, so the degree of crosstalk will vary in time and space.Published spectra show that GRAB ACh3.0 (or EGFP) is essentially not excited at 565 nm, whereas jRGECO1a (or mApple) fluorescence is excited at 470 nm.To quantify the amount of crosstalk under well-controlled conditions, we collected images of purified EGFP and mApple proteins with various combinations of illumination power and exposure time for both the 470-and 565-nm LEDs.In line with the published spectra, we observed that the 565-nm LED practically does not excite purified EGFP while mApple shows fluorescence when excited with the 470-nm LED [Figs.S2(A) and S2(B) in the Supplementary Material].Based on this in vitro calibration with EGFP and mApple, and the typical acquisition parameters we used for in vivo imaging in GRAB ACh3.0 -and jRGECO1aexpressing animals, we determined that the level of crosstalk of GRAB ACh3.0 to fluorescence measured upon 565-nm excitation is <1%, whereas fluorescence from jRGECO1a contributes about 3% to 9% to the fluorescence intensity measured upon 470-nm excitation (Fig. S2 in the Supplementary Material; see Appendix C for a detailed discussion of the crosstalk analysis).To verify that the degree of crosstalk is similar in recordings from awake mice expressing GRAB ACh3.0 and jRGECO1a, we applied a simplified correction approach that does not incorporate dynamic changes of hemoglobin absorption on fluorescence excitation and emission (details in Appendix C).After applying the correction, we used cross-correlation analysis to verify that the degree of crosstalk is within the same range as established under in vitro conditions (Fig. S3 in the Supplementary Material).Here, it should be noted that high correlation between GRAB ACh3.0 and jRGECO1a signals is expected primarily due to their underlying physiological covariance 33 and is only influenced to a small degree by optical crosstalk, as shown in the analysis.However, while the in vitro calibration provides us with the parameters necessary to perform crosstalk correction of in vivo data, we decided not to further use crosstalk correction for the recorded in vivo data as it likely will add additional noise without sufficient benefit-at least under the experimental conditions present in this study (see Appendix C for further discussion).
In addition, it has been reported for mApple-based calcium indicators, such as jRGECO1a, that illumination with blue light leads to photoswitching, which causes a reversible, calciumindependent increase in fluorescence intensity. 46In a control experiment with an awake jRGECO1a-expressing mouse, we did not observe jRGECO1a photoswitching under 470-nm illumination conditions typical for two-color fluorescence imaging in GRAB ACh3.0 -and jRGECO1a-expressing mice.However, when the 470-nm illumination intensity was increased by an order of magnitude compared to typical imaging conditions, we observed photoswitching of jRGECO1a (Fig. S4 in the Supplementary Material).
To demonstrate the capabilities of our new imaging system design under in vivo conditions, we recorded spontaneous and stimulus-induced changes in neuronal calcium, ACh, and hemodynamic signals in awake head-fixed mice expressing jRGECO1a and GRAB ACh3.0 (Fig. 2).Mice were implanted with a curved glass cranial window exposing both hemispheres of the dorsal cortex while leaving the bone in-between the two hemispheres above the superior sagittal sinus intact 48 [Fig.2(a)].Images were registered to the Allen Mouse Common Coordinate Framework 49 to generate time courses from predefined regions of the dorsal cortex.Measurements of spontaneous activity show a high correlation between neuronal calcium and ACh time courses [Fig.2(b) shows right secondary motor cortex], which is consistent with published results. 33We then recorded the response to a sensory stimulus consisting of a 2-s, 3-Hz air puff sequence (20 repetitions, 25-s interstimulus interval) delivered to the right whisker pad [Figs.2(c)-2(e)].Ratio maps show neuronal activation and functional hyperemia in the contralateral barrel cortex [green ROI in calcium map, Fig. 2(c)].The stimulus-averaged time course calculated from the contralateral barrel cortex shows that the imaging system operates at a sufficient signal-to-noise ratio to record the neuronal calcium response to the six individual air puffs comprising one stimulus train [Fig.2(d)].Despite high trial-to-trial variability [Fig.2(e)], which is consistent with published widefield imaging studies in awake mice, 26 responses to individual air puffs are still discernible in individual trial recordings.The time course of the hemodynamic response shows a broad peak and poststimulus undershoot.The hemodynamic correction we apply to the green fluorophore compensates for darkening that occurs when hemoglobin levels increase.Figure S5 in the Supplementary Material compares the trial-averaged ACh time courses with and without correction for hemodynamic artifacts.
With a unit price of >$20,000, the sCMOS camera is the most expensive component of the imaging system.To show how the system performance is affected by camera choice, we repeated the experiment shown in Fig. 2  The aluminum hemisphere in-between the cranial window and the imaging system reflects 525-and 625-nm light and provides diffuse illumination minimizing specular reflections at the curved glass window that interfere with reflectance imaging.Inhomogeneities of signal intensity across the sample surface likely occur due to the curvature of the "crystal skull" glass window [Fig.2(a)].While enabling diffuse illumination, the aluminum hemisphere also acts as a shield preventing the strobing illumination light from reaching the mouse's eyes.For example, fourcolor imaging with an effective acquisition rate of 10 Hz would produce an undesired visual stimulus at 40 Hz.To illustrate the shielding effect, we recorded changes of the animal's pupil diameter during mesoscopic imaging in presence or absence of the aluminum hemisphere (Fig. 3).For this experiment, a behavioral camera recorded the pupil diameter continuously while widefield illumination across the four wavelengths at 40 Hz was switched on and off at 60-s intervals; we recorded 20 off-to-on and on-to-off transitions per test.The images shown in Fig. 3(a) show the pupil before and after turning the illumination on; Video 1 shows a recording of the off-to-on transition with and without hemisphere.Figure 3(b) shows representative 2-min time series covering the off-to-on and on-to-off transitions in presence or absence of the hemisphere.Without the hemisphere, pupil size measurements are largely dominated by the state of widefield illumination.In the presence of the hemisphere, pupil size changes are largely independent from widefield illumination [Fig.3(c)].Then we compared the relative change in pupil size from 1.5 to 0.5 s before versus 2 to 3 s after start of widefield illumination averaged over 20 transitions and observed a significantly larger drop in pupil size without the hemisphere [p < 0.001; t-test; Fig. 3(d)].Pupil tracking is often used as an important variable to define behavioral states 50,51 and it is therefore desirable to combine pupil recordings and mesoscale imaging. 40In our system, the hemisphere eliminates interference of widefield illumination on behavioral recordings and enables high-quality measurements of pupil dynamics during mesoscale imaging.

Discussion and Conclusion
Here, we describe a new widefield microscope design for mesoscopic imaging of dorsal mouse cortex.Several groups describe widefield calcium imaging with a single fluorescent indicator and concurrent hemodynamic imaging, whereas other groups describe two-color fluorescence widefield imaging without hemodynamic imaging; however, we identified one report describing a system for combined two-color fluorescence and hemodynamic imaging using two cameras. 35ompared to previous publications, our system has several unique features: First, the system enables quasi-simultaneous imaging of two spectrally distinct fluorophores (here: EGFP-based GRAB ACh3.0 and mApple-based jRGECO1a) together with quantitative estimates of Δ½HbO and Δ½HbR via two-color absorption imaging on a single detector.Changes in [HbO] and [HbR] are used to correct hemodynamic artifacts of the recorded fluorescence signals.Further, we can use them to study the relationship between cortical hemodynamics and different aspects of brain activity, such as neuronal firing and neuromodulation, thereby gaining further insights into brain state-dependent regulation of cerebral blood flow and metabolism. 52econd, we use a single detector.Using one camera drastically reduces the overall cost to build the system but requires sequential image acquisition for every wavelength (here: four wavelengths).We achieve sufficient SNR with negligible spectral crosstalk at a 10 Hz acquisition rate, which is within the range for widefield imaging studies in awake animals. 23,33inally, the illumination light is contained within an aluminum hemisphere.In addition to providing diffuse illumination for light from 525-to 625-nm LEDs to minimize reflections at the glass window, the hemisphere blocks most of the illumination light from the animal's view.This prevents interference with behavioral measurements of pupil diameter as well as unwanted stimulation of the visual system.Pupil size measurements are often used as a noninvasive, indirect readout of alertness. 50,53We expect our shielded illumination design will enable seamless integration of mesoscale imaging with behavioral experiments based on visual stimulation 27,51 and virtual reality.The overall cost to build our widefield imaging system with the components used is around $45,000, with the camera comprising the bulk of the cost (∼$25,000).Although this is a lower price than for typical commercially available widefield imaging systems, modifications could further reduce the overall cost of the system.Depending on the brightness of the used fluorophores, their expression level, as well as experimental demands, sufficient imaging performance could be achieved with a less expensive camera (see Fig. S6 in the Supplementary Material).Furthermore, we used commercially available LED drivers as well as commercially available LEDs for fluorescence excitation.Using a custom design or already available open-source components with similar performance characteristics, especially regarding their response time to digital triggers, could further reduce the overall cost (see, e.g., Ref. 41).
Our imaging system can easily be modified or expanded to suit the requirements of different experimental studies; changing an LED and the corresponding optical filters will allow combining other fluorescent indicators if they are spectrally separated.This would, for example, allow integration of calcium indicators excited with near-infrared light. 55As the spectral diversity of genetically encoded molecular probes expands, widefield imaging with three spectrally separated fluorescent indicators might soon become a feasible approach.Furthermore, we and others combined mesoscale imaging with extracellular electrophysiology, providing simultaneous access to cortical network activity as well as single-or multiunit activity in targeted cortical and subcortical areas. 21,56,57ur current system design has the following limitations.First, despite aiming to minimize specular reflections and illumination inhomogeneity of 530-and 630-nm LEDs through diffuse illumination, we still observe inhomogeneous levels of absolute image brightness across the field of view.In practice, however, some variations of absolute brightness will not affect the readout of relative (i.e., intensity-normalized) changes in fluorescence or absorption, except for potentially reduced signal-to-noise levels in darker regions of the image.A possible way to further improve illumination homogeneity of both reflectance and fluorescence illumination is the incorporation of an integration device with uniformity correction, such as EUCLID 58 into our imaging system.Further, uneven signal intensity in the medial-to-lateral direction is due to the curvature of the glass window covering the cranial exposure [see Fig. 2(a)].Use of curved glass also causes different cortical depths to match the focal plane of the imaging system, which potentially introduces spatially heterogenous sampling bias.Here, the use of adaptive lens optics such as metalenses 59 could improve uniformity of illumination and could be used to match the focal plane across the field of view.Designing an imaging system with a single camera reduces overall system cost.However, the use of two cameras would allow for higher acquisition rates through spectral multiplexing.Further, the dual-band dichroic mirror and emission filter cut the available spectral bandwidth of the EGFP∕GRAB ACh3.0 emission band (500 to 540 nm) and the mApple/ jRGECO1a excitation band (550 to 570 nm).This potentially enables crosstalk due to off-peak excitation of the red fluorophore (see further discussion in Appendix C).However, with the fluorophores at the expression levels used, we achieve sufficient SNR with negligible spectral crosstalk at a 10-Hz acquisition rate, which is within the range of appropriate temporal resolution for widefield imaging studies in awake animals. 23,33Furthermore, the extension of the imaging system by a second camera requires minimal modifications to the overall layout, allowing for higher acquisition speeds, if desired.Another consideration during the design of our system was the use of lasers instead of LEDs for excitation of the green and red fluorescent proteins.While less expensive, LEDs have a broad emission spectrum that requires clean-up through excitation filters.This limits the emitted light to the usable spectral range [see, e.g., Fig. 1(e), 565-nm LED] leading to an overall reduction in the available power compared to the nominal output (75 mW versus 8.6 W for the 565-nm LED), potentially requiring longer integration times to achieve sufficient SNRs.Here, the use of a laser would enable excitation at higher power levels, if required.A drawback of using lasers instead of LEDs would be a substantial increase in the cost of the system.
In summary, given the versatility of our system and the comparatively low cost in contrast to commercial intravital imaging setups, we expect our widefield imaging system design to become a useful platform for neurovascular studies in the neuroscience and neurophotonics communities.
4 Appendix A: Methods

Experimental Animals
All procedures were conducted with approval from Boston University Institutional Animal Care and Use Committee and in accordance with the National Institutes of Health Guide for the Care and Use of Laboratory Animals.Standard rodent chow and water were provided ad libitum.Mice were housed in a 12-h light cycle with lights turning on at 7:30 am.We used three female mice of the Thy1-jRGECO1a GP8.20 16,46 line; cortex-wide GRAB ACh3.0 expression was induced by injection of 1.5 μL AAV9-hSyn-GRAB ACh3.0 (WZ Biosciences, 3.06 × 10 13 GC∕mL) into each transverse sinus (3 μL per animal) of 1-day old neonates according to procedures described by Lohani et al. 33 Headpost implantation and craniotomy in adult (8-to 16-week old) animals follows procedures described by Kilic et al. 48Briefly, a custom-designed headpost machined from titanium was attached to the cranium and a modified crystal skull curved glass window was used to replace the dorsal cranium on both hemispheres. 10Prior to surgery, the original curved glass (12 mm width, labmaker.org)was cut in half to obtain separate glass pieces for each hemisphere.Following the implantation procedure, a silicone plug was placed on top of the glass window and a 3D-printed cap was fixed to the headpost to prevent heat loss.Mice were allowed to recover for at least 1 week following surgery before behavioral training commenced.
To adapt mice to be comfortable under head fixation, animals were handled until they remained calm in the gloved hand of the investigator.Mice were then head-fixed for increasing periods while they received a reward of sweetened condensed milk.Data acquisition was started when animals tolerated 1 h of head fixation.

Imaging System Optical Components
The imaging system is assembled from standard optomechanical components compatible with 2inch lens tubes.The imaging system is placed on an elevated platform (61 × 61 cm, Thorlabs, MB2424) on a vibration-isolating table.
An objective lens (Olympus, ZDM-1-MVX063) focuses the fluorescence excitation light onto the brain for epi-illumination and projects the emitted light toward the camera.The objective lens is surrounded by the RIR assembly, which contains six equally spaced 525-nm LEDs (CREE LED, XBDGRN-00-0000-000000C02) and six equally spaced 625-nm LEDs (CREE LED, XBDRED-00-0000-000000701) for reflectance imaging.Bandpass filters were cut to 9.4-mm diameter and placed in front of the reflectance LEDs.These filters are 532/10 nm bandpass filters (Edmund Optics, 65094) for the 525-nm LEDs and 632∕10 nm bandpass filters (Edmund Optics, 65105) for the 625-nm LEDs.Six Peltier elements (Ferro Tec, 81044) controlled by a feedbackcontrolled temperature controller (Thorlabs, TED200C) actively cool the LEDs in the RIR to maintain the LEDs at room temperature during data acquisition.Temperature is measured on the LED ring with a 10-kΩ thermistor (TDK, NTCG163JF103FT1).Heat from the Peltier elements is transferred to a custom heat sink made of aluminum.The RIR allows fixation of the custom aluminum hemisphere (bead-blasted) with a 15-mm opening located directly above the cranial window.
To perform imaging tests with a different camera model (ac1920-150uc, Basler with Python 2000 CMOS chip from ON semiconductor in global shutter mode), a 45-deg mirror (Edmund Optics, 83-538) behind the detection filter was used to divert the light path to a second camera arm with alternative tube lens (ThorLabs, TTL180-A).The camera was controlled through the manufacturer's software (Pylon Viewer); 4 × 4 binning (sum) with analog gain disabled was used with exposure times set to 6 ms for all frames.

Control of Imaging Acquisition
We use a computer-controlled DAQ device (National Instruments, PCIe-6363) to operate the imaging system.The DAQ device is connected to two I/O connector boards (National Instruments, CB-68LPR) inside a control interface box.The control interface relays a DAQgenerated exposure trigger to the camera.An IC with an OR gate (Texas Instruments, SN74AHC32N) within the control interface receives the "Fire ALL" trigger from the camera as an input.The Fire ALL trigger is positive when the exposure state overlaps for all rows in the rolling shutter mode.The output of the OR gate is the common input to an IC with four AND gates (Texas Instruments, SN74AC08N).The second-wavelength-specific-input to each individual AND gate is the LED selection trigger generated by the DAQ system.Each output of the AND gate connects to one LED driver (Thorlabs, DC2200) controlling one type of LED.During normal operation of the imaging system, the "Fire ALL" trigger controls when the LED turns on.The other input to the OR gate is a camera bypass trigger that allows the LEDs to be triggered while the camera is off.We use a custom MATLAB script to control the DAQ system.We used image capture software (Andor, SOLIS) to control camera and image acquisition.
The animal was positioned under 525-nm illumination and the platform carrying the animal was adjusted to bring the brain surface into focus.Then the aluminum dome is lowered to minimize the distance between its opening and the brain surface.Illumination intensities (i.e., LED currents) were adjusted to bring image intensities into the dynamic range of the camera at 16-bit image depth.The acquisition window was cropped to the size of the field of view; image height (number of rows) was a variable entered into the MATLAB script controlling image acquisition: the length of the period when all rows are exposed ("FIRE all") depends on image height and requires adjustment of the camera exposure trigger generated by the DAQ system.Typical image sizes were 460 × 360 pixels with a 9.2 × 7.2 mm field of view and 20-μm resolution at 4 × 4 binning.Based on selected exposure times and frame rate, the MATLAB script calculates the total number of frames per acquisition run; this number is entered into the image capture software.The camera is set to "External Exposure" mode and image acquisition is started.Then execution of the MATLAB script is started causing the DAQ interface to send LED selection and camera triggers to the control interface.

Sensory Stimulation and Behavioral Recordings
The MATLAB acquisition script generates additional triggers for synchronization with behavioral apparatuses, such as to deliver stimuli by puffing compressed air to the animals' whiskers.Stimuli are defined by duration, frequency, repetitions, interstimulus interval for individual stimuli and sequences of prestimulus baseline, stimulus, and poststimulus observation period.Based on these parameters, the MATLAB script calculates the required imaging time and the number of images.For air puff stimuli, the 5-V trigger is sent to a picopump (WPI, PV830) connected to a 2-mm diameter plastic tube that terminates in a glass capillary to deliver pressurized air to the mouse whiskers.
The animal's face and whiskers are recorded during imaging with a CMOS camera (Basler, acA1920-150uc) through a variable zoom lens (Edmund Optics, 67715).A 940-nm LED (Thorlabs, M940L3) continuously illuminates the animal, and a 920-nm longpass filter prevents light from the microscope from contaminating the camera recording.The mesoscope control system triggers the CCD camera to acquire a single frame during each acquisition cycle.
An accelerometer (Analog Devices, ADXL335) is placed underneath the mouse bed to record animal's motion.The accelerometer generates analog signals for motion in X, Y, and Z direction, which are recorded by a separate DAQ system (NI, USB-6363) receiving a synchronization signal from the mesoscope control system.

Image Processing
The camera control software (Solis, Andor) stores the acquired images as raw data (".dat")together with metadata files (".ini" and ".sifx").For all subsequent import and processing steps, we use MATLAB.Raw image data are loaded into MATLAB, bypassing conversion into image (.tif) files by the camera software and enabling multistream data import from high-performance solid-state drives with MATLAB's multiprocessor environment.For all channels, a baseline image representing the temporal average of all images for that channel is generated.For fluorescence channels, ΔF∕F images are computed on a pixel-by-pixel basis by dividing each image by the baseline image.After calculating ΔF∕F, linear trends are removed from the time course of each pixel using the MATLAB function "detrend."The modified Beer-Lambert law is used to calculate changes in HbO and HbR concentrations from 525 and 625 nm reflectance (see implementation for MATLAB's multiprocessor environment in Appendix B).Estimated Δ½HbO and Δ½HbR are used to correct green fluorescence measurements for artificial intensity changes caused by dynamic changes in excitation and emission light absorption by hemoglobin.The correction method has been described by Ma et al. 4 and its implementation is further detailed in Appendix B. In this correction approach, the measured fluorescence is equal to the artifact-free change in fluorescence multiplied by a correction factor containing the pathlength in tissue and the tissue absorption coefficient.Pathlengths have been previously estimated based on Monte-Carlo simulations. 4Since hemoglobin accounts for most of the light absorption in brain tissue, light absorption depends largely on Δ½HbO and Δ½HbR and the known HbO and HbR extinction coefficients at the respective wavelengths of fluorescence excitation and emission.We calculate the light absorption coefficient for every pixel at every time point based on our estimations of Δ½HbO and Δ½HbR and apply it to correct for hemodynamic darkening of ΔF∕F in the green (EGFP∕GRAB ACh3.0 ) channel.
Measurements with several consecutive trials are interpolated at 100-ms intervals using a time vector relative to the stimulus onset.Trial-averaged ratio maps are calculated by averaging individual time courses across trials.
Behavioral video recordings are stored as image sequences and imported into MATLAB.To calculate pupil size, an elliptical region is drawn defining the eye and an intensity threshold is manually defined.The pupil size is calculated as the number of pixels below the threshold.The resulting time series is normalized to its maximal and minimal pupil size within the recording.

In Vitro Crosstalk Calibration
Solutions of purified 0.03 mg∕mL EGFP (Abcam, ab84191) and 0.4 mg∕mL mApple (kindly provided by Ahmed Abdelfattah) in phosphate-buffered saline (PBS) were prepared prior to imaging.30-μL samples of EGFP, mApple, or PBS-only control solution were each injected into one channel of separate 6-channel slides (ibidi, μ-Slide VI-Flat).The center of the fluorophore-filled channel was focused under the widefield imaging system before starting image acquisition for each slide.Illumination during an imaging sequence alternated between the 470-and 565-nm LEDs while stepping through a series of camera exposure times and increasing LED currents.Three trials were conducted for each sample, with the position of the slide under the objective being shifted between trials to account for intensity variations.To calculate the fluorescence intensity, pixel values were averaged over a region of interest covering the center of the fluorophore-containing channel.Then we took the mean of the three trials per condition and subtracted the fluorescence intensity of the PBS control.Further analysis of crosstalk is detailed in Appendix C.

Photoswitching of jRGECO1a
A variant of the MATLAB acquisition script was used to create two types of acquisition cycles: One cycle includes 470-nm illumination but disables image acquisition during 470-nm illumination (no camera exposure trigger and the FIRE ALL trigger is bypassed), whereas the other cycle does not include 470-nm illumination and has normal image acquisition.The combination of both cycles enables intermittent 470-nm illumination cycles at variable frequencies.Image analysis is performed on baseline-normalized time courses averaged from all pixels across the brain surface.) can be estimated using a system of linear equations (here, we use intensity measurements at two wavelengths (525 and 625 nm) and two unknowns ([HbR] and [HbO]).For further discussion, see Ma et al. (2016). 18he molar extinction coefficients ε of Hb and HbO at 525 and 625 nm are known variables: E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 0 1 ; 1 1 7 ; 6 4 3 E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 0 2 ; 1 1 7 ; 6 0 7 We estimate the optical density in every pixel for every time point t from the measured intensities I t 525 and I t 625 as follows: E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 0 3 ; 1 1 7 ; 5 7 5 E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 0 4 ; 1 1 7 ; 5 1 9 The pathlengths X 525 and X 625 are known variables; the "baseline" intensities I 0 525 and I 0 625 are estimated using the temporal average of the intensities at the given pixel.
When we solve Eqs. ( 1) and ( 2) explicitly, combine them with Eqs. ( 3) and ( 4), and combine all variables that are constant over space and/or time, we enable an efficient (i.e., fast) approach for estimating [Hb] and [HbO].
To estimate Δ½HbO, we first solve Eqs. ( 1) and ( 2) for Δ½Hb: Then we solve Eq. (5c) for Δ½HbO: Then we combine all constant variables in Eq. (5g) to define C HbO 525 and C HbO 625 : E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 0 5 h ; 1 1 7 ; 1 3 2 Now, we insert the definition of ΔOD 525 and ΔOD 625 from Eqs. ( 3) and ( 4): ; t e m p : i n t r a l i n k -; e 0 0 6 a ; 1 1 7 ; 9 4 and combine again all variables that are constant over time:  3) and ( 4): E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 0 9 a ; 1 1 4 ; 1 0 2 Then, we combine again all variables that are constant over time: Note that C HbÃ 525 and C HbÃ 625 are invariant in time and space, whereas A Hb 0 is only time-invariantits value changes for every pixel.
The implementation in MATLAB is as follows: ( We performed fluorescence measurements of purified EGFP and mApple that were excited at 475 nm and 565 nm using series of increasing exposure times and optical powers.For further evaluation, we calculated the radiant energy E for each measurement as the product of illumination power, estimated from the LED current using a separately measured power calibration curve, and exposure time. Then we defined the illumination ratio R for EGFP and mApple as radiant energy at off-peak illumination divided by radiant energy at on-peak illumination combining different illumination conditions for on-and off-peak excitation into a single variable: E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 1 1 ; 1 1 7 ; 1 8 0 (11)   and E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 1 2 ; 1 1 7 ; 1 4 4 Note that every illumination ratio can be achieved by different combinations of optical powers and exposure times.This means that pairs of fluorescence measurements are performed with different combinations of optical power and exposure time for off-and on-peak excitation but match the chosen illumination ratio.The corresponding individual on-and off-peak fluorescence intensities have different values.We define their ratio, measured at a certain LED power P 1 and exposure time t 1 at off-peak excitation relative to the fluorescence intensity at a certain LED power P 2 and exposure time t 2 at on-peak excitation for each fluorophore as crosstalk ratio C for EGFP and mApple as E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 1 3 ; We again substitute LED power and exposure time by radiant energy E: To estimate the crosstalk ratio C, we chose 20 pairs of measurements matching a given illumination ratio R. The pairs of fluorescence intensities from off-peak excitation versus on-peak excitation, when plotted against each other, follow a linear relationship [Figs.S2(C) and S2(D) in the Supplementary Material], which is expected given the underlying linear relationship between fluorescence intensity and radiant energy [shown in Figs.S2(A) and S2(B) in the Supplementary Material].Accordingly, we can estimate the crosstalk ratio C from the slope of the linear fit between fluorescence from off-and on-peak excitation for the 20 pairs of values for each given illumination ratio [Figs.S2(C) and S2(D) in the Supplementary Material].Finally, for each fluorophore, we plot the crosstalk ratio C versus illumination ratio R for EGFP and mApple (i.e., C EGFP versus E 565 nm ∕E 470 nm and C mApple versus E 470 nm ∕E 565 nm ; Fig. S2(E) in the Supplementary Material).In the resulting graph, a small value of R indicates weak off-peak excitation compared to on-peak excitation, and the resulting crosstalk ratio C is small.For increasing values of R, off-peak excitation becomes increasingly stronger compared to on-peak excitation leading to a higher crosstalk ratio C. Given that EGFP shows only very weak off-peak excitation at 565 nm, its crosstalk ratio C increases only modestly with R. In contrast, as mApple is more efficiently excited off-peak at 470 nm, the crosstalk ratio C rises faster with increasing R.

6.2
In-Vivo Correction of jRGECO1a Crosstalk Upon 470-nm Excitation in GRAB ACh3.0 and jRGECO1a-Expressing Animals In experiments with jRGECO1a-and GRAB ACh3.0 -expressing mice (Fig. 2), the used LED powers and exposure times for fluorescence excitation at 470 and 565 nm define a range of illumination ratios.Based on the in vitro calibration with EGFP and mApple, we estimated that for EGFP-derived GRAB ACh3.0 , <1% of its fluorescence present at 470-nm excitation contributes to fluorescence measured at 565-nm excitation, whereas for mApple-derived jRGECO1a, 3% to 9% of its fluorescence present at 565-nm excitation contributes to fluorescence measured at 470-nm excitation [highlighted areas in Fig. S2(E) in the Supplementary Material].
Calibration of the crosstalk ratio C enables us to estimate the level of fluorescence present at 470-nm excitation due to off-peak excitation of jRGECO1a for a given combination of illumination power and exposure times, i.e., radiant energies E 1 and E 2 : We assume here that fluorescence measured at 565 nm is solely due to jRGECO1a fluorescence with negligible contribution from GRAB ACh3.0 , given the minor (<1%) contribution of EGFP fluorescence at 565-nm excitation.This enables us to define the following correction term to remove jRGECO1a crosstalk from fluorescence measured at 470-nm excitation: E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 1 8 ; 1 To apply this correction approach based on in vitro determined crosstalk ratios C to in vivo data, we make the following assumptions: (a) fluorescence acquired at 565 nm is solely caused by jRGECO1a and (b) dynamic changes in jRGECO1a fluorescence between consecutive acquisition of images at 470-and 565-nm excitation are negligible.Further, we do not consider wavelength-dependent absorption of light for fluorescence excitation and emission through hemoglobin.Incorporating [HbO] and [HbR] estimates into the correction would result in combining measurements from all four wavelengths to correct for crosstalk in fluorescence measurements at 470-nm excitation, potentially confounding this measurement, i.e., GRAB ACh3.0 fluorescence, by the sum of measurement noise from all measured channels.
To verify that the degree of in vivo crosstalk matches our in vitro estimation, following Eq.( 19), we assume that we can estimate the "true" GRAB ACh3.0 signal by subtracting a proportion C of the 565-nm fluorescence image intensity from the 470-nm fluorescence image on a pixel-by-pixel and frame-by-frame basis.Here we subtracted three different crosstalk ratios (we tested C ¼ 0, 3%, 5%, and 10%) and estimated the cross correlation of the resulting, corrected F 470 timeseries with the F 565 timeseries (Fig. S3 in the Supplementary Material).For C of 5% and 10%, we observe clear, negative zero-lag cross-correlation indicative for overcorrection.This result suggests that the degree of crosstalk established under in vitro conditions is a reasonable estimate of the crosstalk we observe in vivo.Based on this evaluation, we suggest that under the illumination conditions for in vivo imaging used here, the benefits of spectral crosstalk correction are minimal.Therefore, we did not perform crosstalk correction for the data shown in Fig. 2.
However, it is possible that levels of crosstalk increase, particularly when fluorescence from on-peak excitation of one fluorophore approaches levels of the other fluorophore due to off-peak excitation at the same wavelength.First, this could happen due to low expression of one fluorescent indicator versus the other.Here, using AAV delivery into newborn pups, we achieved relatively high levels of GRAB ACh3.0 expression, diminishing the potential effect of jRGECO1a crosstalk on GRAB ACh3.0 measurements.Second, this situation could also arise when one fluorescent indicator has low levels of basal fluorescence in virtual absence of its analyte, i.e., ACh for GRAB ACh3.0 in this study.When ACh levels are low (low GRAB ACh3.0 fluorescence) and neuronal calcium levels are high (high jRGECO1a fluorescence), the measurement at 470-nm excitation could potentially show higher contribution from jRGECO1a fluorescence compared to periods where ACh levels are high.To further address and quantify this effect, however, is very difficult with the used animals under the given imaging conditions as jRGECO1a and GRAB ACh3.0 signals show high correlation due to their underlying physiological covariance.If further examination of jRGECO1a crosstalk into the green fluorescence channel is required, animals expressing jRGECO1a and EGFP, or a mutated version of a green-fluorescent indicator that does not show activity or ligand-dependent changes in fluorescence could be used.
review, and editing.Dora Balog performed software, review, and editing.Bernhard Zimmermann, Kıvılcım Kılıç, Joel D. Herbert, Grace Chabbott, Bradley C. Rauscher, and John X. Jiang performed resources, review, and editing.Sreekanth Kura and Harrison P. Fisher performed software, review, and editing.Sava Sakadzic and David A. Boas performed conceptualization, review, and editing.Anna Devor performed conceptualization, funding acquisition, supervision, writing of the original draft, review, and editing.Ichun Anderson Chen performed conceptualization, methodology, supervision, review, and editing.Martin Thunemann performed conceptualization, software, supervision, visualization, writing of the original draft, review, and editing.

Fig. 1
Fig. 1 Schematics of the wide field imaging system and fluorescence and reflectance spectra.(a) A single sCMOS camera is used to measure emission of green and red fluorophores and two hemoglobin absorption wavelengths.A 565-nm LED excites mApple-based jRGECO1a and a 470-nm LED excites EGFP-based GRAB ACh3.0 ; both LEDs are filtered through excitation filters.The main dichroic mirror reflects excitation light toward a 45-deg-mirror and then a 0.63× objective lens to focus the excitation light onto the sample.Below the objective lens is the reflectance illumination ring (RIR) further detailed in panel (b).An aluminum hemisphere between the RIR and the sample surface internally reflects the 525-and 625-nm LEDs for diffuse illumination and prevents specular reflections and illumination light from reaching the animal's eyes.Fluorescence emission and reflectance pass the main dichroic and the detection filter before being focused by a tube lens to the camera.(b) Cross section of the RIR setup.Each of the six 525-and 625-nm LEDs are arranged on a ring-shaped board at equal distances.Peltier elements actively cool the LEDs to prevent temperature changes resulting in unstable illumination intensities.10-nm wide excitation filters are mounted in front of each LED.(c) Timing diagram to trigger LEDs and camera acquisition: MATLAB controls the DAQ system via digital triggers controlling camera and LED drivers.A different LED illuminates the sample during each frame for near-simultaneous acquisition.The effective frame rate is determined by the length of the four-frame cycle (typically 100 ms).Grey lines next to the camera Fire ALL trigger represent exposure of pixel rows in rolling shutter mode.The Fire ALL trigger is sent when exposure overlaps for all rows.Each LED driver (bottom left) is controlled by the output of a separate AND gate.The Fire ALL trigger is the common input to every AND gate, whereas the individual LED selection triggers are the second input for each AND gate.(d).Top: light-blue and dark-blue shaded areas represent the 470-nm LED spectrum before and after passing the excitation filter (BP466/40, blue dotted line).Bottom: excitation (solid blue line) and emission (solid green line) spectrum of EGFP-derived GRAB ACh3.0 . 17(e) Top: light-orange and dark-orange areas represent the 565-nm LED spectrum before and after passing the excitation filter (BP560/14, orange dotted line).Bottom: excitation (solid yellow line) and emission (solid red line) spectrum of mApple-derived jRGECO1a.16(f) Top: LED spectra are confined by excitation filters (525 nm in light and dark green; 625 nm in light and dark red before and after passing through the respective excitation filters).Bottom: absorption spectra of oxyhemoglobin (HbO, red line) and deoxyhemoglobin (HbR, blue line) from Prahl.47 (g) Spectra of the main dichroic (brown) and detection filter (pink).Light must pass through both filters to reach the camera.

Fig. 2
Fig. 2 Spontaneous and stimulus-evoked activity in the somatosensory cortex.(a) Average of the first ten images collected during a 10-min acquisition period for each illumination wavelength.The two hemispheres were manually masked around the cranial window.A green polygon was drawn around the right secondary motor cortex, which was determined by registration to the Allen Atlas. 49(b) Spontaneous activity in the right secondary motor cortex [green polygon in panel (a)] for calcium, acetylcholine (ACh), and oxygenated, deoxygenated, and total hemoglobin (HbO, HbR, and HbT).(c) Ratio maps averaged across 20 trials showing the response to a 2-s, 3-Hz train of air puffs to the right whisker pad.Maps show data from 1.7 s after stimulus onset except for the calcium map, which shows data from 0.1 s after stimulus onset.(d) Average time course of the stimulus-evoked response in the contralateral barrel cortex [green region in panel (c)].The pink dashed line indicates when the calcium ratio map in panel (c) is shown, whereas the black dashed line indicates when ACh, HbT, and HbO ratio maps in panel (c) are shown.The grayshaded area indicates the duration of the stimulus.(e) Response of the contralateral barrel cortex to individual stimulus trains.The trials are sorted according to the magnitude of the calcium response.
with a CMOS camera costing around $750 (Fig. S6 in the Supplementary Material).Time courses of signals averaged from spatially connected pixels in a region of interest show stimulus-evoked and spontaneous changes in neuronal calcium, ACh, and hemodynamic signals [Figs.S6(B), S6(D), and S6(E) in the Supplementary Material].On the other hand, we observed spatially structured noise across the field of view [compare Fig. 2(c) with Fig. S6(C) in the Supplementary Material] that could impact pixel-based analysis methods.

Fig. 3
Fig. 3 Aluminum hemisphere shields illumination and reduces pupil constriction during mesoscale imaging.(a) Images of the mouse face 0.5 s before and 5 s after widefield imaging acquisition begins.Green circles highlight the pupils.(b) Example 2-min pupil diameter time courses showing the off-to-on and on-to-off transition with (top) and without (bottom) hemisphere.The gray-shaded area indicates the period when widefield acquisition was performed.(c) Average pupil diameter changes upon beginning of widefield imaging (20 trials; traces show average ± standard deviation).Pupil diameter 1.5 to 0.5 s before imaging was defined as baseline for each trial.The green-dashed line indicates the start of widefield imaging, the black dashed lines indicate time points of the images shown in panel (a).(d) Comparison of the ratio of pupil size 2 to 3 s after start of imaging to pupil size 1.5 to 0.5 s before imaging.Significance was tested with a student's t test (Video 1, MP4, 5.93 MB [URL: https://doi.org/10.1117/1.NPh.11.3.034310.s1]).

5
Appendix B: Optimization of Δ[HbO] and Δ[HbR] Estimation Process for MATLAB Using an Explicit Solution of HbR and HbO with Measured Reflectance at 525 and 625 nm Concentration changes of deoxy-and oxyhemoglobin ([HbR] and [HbO] E T ; t e m p : i n t r a l i n k -; e 0 T A R G E T ; t e m p : i n t r a l i n k -; e 0 .(5a) and (5b):E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 0 5 c ; 1

e m p : i n t r a l i n
T A R G E T ; t e m p : i n t r a l i n k -; e 0

E
Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 T A R G E T ; t e m p : i n t r a l i n k -; e 0 combine all constant variables into three parameters A HbO 0 , C HbOÃ 525 , and C HbOÃ 625 :E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 T A R G E T ; t e m p : i n t r a l i n k -; e 0 T A R G E T ; t e m p : i n t r a l i n k -; e 0 T A R G E T ; t e m p : i n t r a l i n k -; e 0 0 7 d ; 1 time and space, whereas A HbO 0 is only timeinvariant-its value changes for every pixel.In the same way, we can estimate[Hb].First, we solve Eqs.(1) and (2) for [HbO]:E Q -T A R G E T; t e m p : i n t r a l i n k -; e 0 T A R G E T ; t e m p : i n t r a l i n k -; e 0 .(8a) and (8b):E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 0 8 c ; 1We solve Eq. (8c) for Δ½Hb:E Q -T A R G E T ; te m p : i n t r a l i n k -; e 0 0 8 d ; 1 T A R G E T ; t e m p : i n t r a l i n k -; e 0 We combine all constant variables in Eq. (8g) to define C Hb 525 and C Hb 625 E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 0 8 h ; 1 1 4 ; 1 3 8 Δ½Hb ¼ C Hb 625 ΔOD 625 − C Hb 525 ΔOD 525 (8h) and insert the definition of ΔOD 525 and ΔOD 625 from Eqs. (

E
Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 T A R G E T ; t e m p : i n t r a l i n k -; e 0 combine all constant variables into three parameters A Hb 0 , C HbÃ 525 , and C HbÃ 625 :E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 T A R G E T ; t e m p : i n t r a l i n k -; e 0 T A R G E T ; t e m p : i n t r a l i n k -; e 0 T A R G E T ; t e m p : i n t r a l i n k -; e 0 1 0 d ; 1